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Abstract 

Quintessence, a time-varying energy component that may account for the accelerated expansion 
of the universe, can be characterized by its equation of state and sound speed. In this paper, we 
show that if the quintessence density is at least one percent of the critical density at the surface of 
last scattering the cosmic microwave background anisotropy can distinguish between models whose 
sound speed is near the speed of light versus near zero, which could be useful in distinguishing 
competing candidates for dark energy. 
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Recent evidence suggests that most of the energy density of the universe consists of a 
dark energy component with negative pressure. Two candidates are a cosmological constant 
(or vacuum density) and quintessence, [|^ a time-varying energy component with negative 
pressure. A common example of quintessence is a scalar held slowly rolling down its po¬ 
tential. The scalar held may be regarded as real, or simply as a device for modeling more 
general cosmic huids with negative pressure. A cosmological constant and quintessence can, 
in principle, be distinguished by the equation of state, tc, the ratio of the pressure p to the 
energy density p. A cosmological constant always has w = —1 whereas a scalar held gen¬ 
erally has w diherent from unity and time-varying. Through measurements of supernovae, 
large-scale structure and the cosmic microwave background (cmb) anisotropy, the equation 
of state may be determined accurately enough in the next few years to hnd out whether w 
is actually diherent from —1 or not. 

For a scalar held with a negative ehective pressure and slowly varying energy density, 
we can further dehne two cases. The hrst consists of standard models of quintessence |l|, 
in which the scalar held 0 has a canonical kinetic term, X = A second category 

consists of models in which the kinetic energy is not canonical. Prominent examples are 
A:-essence models, which are designed to address the issue of why cosmic acceleration has 
recently begun. The equation of state in A:-essence models is positive until the onset of 
matter-domination triggers a change to negative pressure. The models rely on dynamical 
attractor behavior due to the non-canonical kinetic energy density. The Lagrangian density 
is a rather generic function of X\ C^{X^ cj)) = p[X)/(p. 

One diherence between standard quintessence and k-essence models is the time-evolution 
of the equation of state. The equation of state for a k-essence component approaches — 1 
soon after the onset of matter-domination and then increases towards a less negative value 
in the present epoch as the component begins to dominate the energy density. In stan¬ 
dard quintessence models, the equation of state is generically monotonically decreasing and 
approaching —1 today. This feature was discussed in detail in Ref. 3. 

In this paper, we focus on a second physical property - the speed of sound - which also 
distinguishes standard quintessence from /c-essence and, more generally, from other cosmic 
fluids described by a non-canonical kinetic energy density. Let us consider a scalar held 0 
with a general Lagrangian density Crt>{X,(j)). The stress-energy tensor in this case can be 
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expressed in the form of an ideal fluid 


= {p + p)u>^u^ 


( 1 ) 


where the pressure is p = C^, the energy density is p = and ,X denotes the 

partial derivative with respect to X. From this, we hnd 


w = 


( 2 ) 


The effective speed of sound entering the equations which describe the evolution of small 
perturbations is 

^4>,x 


2 P,x 

Cs = - 


( 3 ) 


P,X ^<I>,X + 2jC^xxX 
Thus, for the standard quintessence models with canonical scalar helds, C = X — V{(j)), the 
equation of state is 

X - 1/ 

w = 


X + V 

and the speed of sound is always equal to one: = 1. However, if C^^xx 7^ 0, then 

7^ 1. (In fact, even > 1 is possible. Physically, this means that perturbations of 
the background scalar held can travel faster than light as measured in the preferred frame 
where the background held is homogeneous. For a time dependent background held, this is 
not a Lorentz invariant state. There is no violation of causality. The underlying theory is 
manifestly Lorentz invariant and it is not possible to transmit information faster than light 
along arbitrary space-like directions or create closed time-like curves.) 

This paper investigates how the variable speed of sound inhuences the huctuations of 
the cosmic microwave background compared to the case of standard quintessence where 
Cg = 1. In general, the ehect is small, but we show that it is detectable in cases like k- 
essence models in which the speed of sound is nearly zero during most of the period between 
last scattering and the present epoch. This behavior produces the greatest diherence from 
standard quintessence. P] We hrst compare models with exactly the same equation of state 
as a function of redshift, w{z), but different sound speed. We hnd that models with near¬ 
zero sound speed today (such as /^-essence models) are distinguishable from models with 
Cg = 1 based on measurements of the CMB power spectrum, provided the quintessence 
energy density is greater than a few percent of the critical density at last scattering. The 
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FIG. 1: The equation of state for the fiducial model (solid line) and the closest spline fit (dashed 
line) quintessence model. The spline fit is chosen so that the CMB spectrum matches the hducial 
model to within the cosmic variance limit. However, the spline itself is smoother than the actual 
w{a). As discussed in the text, there is a large degeneracy in the spline parameters, so the spline 
need not mimic the fiducial equation of state very closely to obtain a good match to the CMB 
spectrum. 

density requirement, which is satished by typical A:-essence models, for example, is needed 
so that the sound speed has a measurable effect on the acoustic oscillation peaks of the 
CMB, which are sensitive to conditions at the last scattering surface. Similar results can be 
obtained for more general forms of dark energy. [^, We then consider whether the effect can 
be mimicked by varying other cosmic parameters or by introducing a time-varying equation- 
of-state. To perform the studies, we introduce a spline technique that is useful in exploring 
models with time-varying w. Our conclusion is that the sound speed effect is distinguishable 
from all other standard parameter effects. Hence, the CMB can provide a useful constraint 
on the sound speed of dark energy. 

In the case of /^-essence the Lagrangian generally has the form = p{X)/cfP‘, where 
p{X) is generally some function with p^xx 7^ 0. We have chosen a specihc form for p for 
our hducial model. The equation of state and sound speed for this model are shown in 
Figs m and H respectively; they are expressed as functions of the scale factor a (a = 1 
today) by integrating the equations of motion. To completely specify the model, we hx the 
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FIG. 2; The speed of sound as a function of the scale factor o (otoday = !)• Note that at <C 1 
near the last scattering surface and at the present epoch. 

cosmological parameters today to typical values: fib = 0.05, flcDM = 0.3, Qq = 0.65 and 
ho = 0.5 (where ho is the Hubble parameter in units of 100 kmsec“^ Mpc“^). The energy 
density as a fraction of the critical density, Hg (A:-essence), is shown in Fig. |[ While c2 > 1 
for large redshifts, this is not an important feature of the model: fig is reasonably small 
whenever > 1, so the value of Cg at these times has negligible effect on the CMB. (We have 
verihed this by rerunning the calculation after artihcially truncating the speed of sound at 
Cs < 1 and comparing the CMB power spectra. Also, we should note that, with a slightly 
different parameterization, we can obtain a model in which Cs < 1 which at early times and 
which has the same behavior at late times.) The most important feature of the hducial 
model is that 1 whenever k-essence contributes signihcantly to the energy density of 

the universe. 

The CMB power spectra for our models are computed by modifying the standard CMB- 
FAST code. i. i Simply put, the effect of the speed of sound on the CMB perturbation 
equations is such that for <C 1, A:-essence will collapse via gravitational instability into 
cold dark matter (cdm) potentials, whereas in the superluminal case, > 1, the growth of 
density perturbations is suppressed. For those familiar with the code, the modihcations are 
straightforward. The perturbed line element is 

ds^ = a‘^{ri)[dri‘^ — {6ij + hij)dx^dx^] (4) 
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FIG. 3: Ofc-essence and Oq as a function of z for the ^-essence (solid line) and best-fit spline models 
(dashed line), respectively. Note that Qq falls off to zero at large redshifts, whereas Qk-essence 
approaches a finite value; w approaches 1/3 for this model at large redshifts. 


where <5^ is the unperturbed spatial metric, and hij is the metric perturbation. We shall use 
h to represent the trace of the spatial metric perturbation. The effect we are examining is 
due to the perturbations to the /c-essence stress-energy in the synchronous gauge for a mode 
with wavenumber k 


dp = -2p— - (p + p)—c, 

4> y 

dp = -2p— - (p + p) — 

0 y 

e = ( 5 ) 

where y = l/\fX and 9 is the divergence of the fluid velocity. The density contrast, d = dp/ p, 
obeys the equation 
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-(1 + w ) 


e + hi 


-m 



( 6 ) 


where the derivative is with respect to conformal time and = a/a. Notice that the 
quantity dp/dp appears, which is not to be confused with c^. In fact 


dp 


^ [37f(l + w){cl -«;)+«;] 


( 7 ) 
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which leads to a simplihed evolution equation for the velocity gradient 


d = {3cl - l)n9 + c^k'^S/il + w). 


( 8 ) 


The kinetic quintessence is distinct from scalar held quintessence, for which = 1 in equa¬ 
tions ( 0 ) and (^. Linearized perturbations in fc-essence can propagate non-relativistically, 
with Cg -C 1. We can see in equation that a small sound speed will cause the velocity gra¬ 
dient to decay; with the conventional gauge choice that Ocdm = 0, the inhomogeneities in the 
fc-essence will describe a huid which is comoving with the cold dark matter. From equation 
(i. we see that the second term on the RHS will be negligible even on scales approaching 
the horizon. The overall effect is that the pressure huctuations 6p are too weak to prevent 
fc-essence collapse via gravitational instability into the CDM gravitational potentials. 

The CMBFAST code takes w{a) and Cs(a) as inputs, so it is possible to manually adjust 
these functions to have any values (including, of course, = 1). Once we have computed 
the CMB anisotropy for two models, they can be compared by computing their likelihood 
difference, the probability that they could be confused due to the cosmic variance in local 
measurements of the CMB. Given models A and S, the negative log-likelihood, — logL, is 
derived in [Q: 


- logL = + I) (l - + log 


j ClB CiB ' 

The condition — log L > 6 corresponds to distinguishability at the 3a level or greater. The 
relative normalization of the spectra is chosen so as to minimize the likelihood difference, 
including I up to 1500. 

The power spectrum for our hducial model is given by the solid line in Fig. We can see 
the effect of the unusual speed of sound on this model by computing a new spectrum, which 
has the same equation of state and cosmological parameters, but whose speed of sound has 
been set equal to 1 for all a. This corresponds to a quintessence held - with canonical kinetic 
energy - rolling down a potential. The power spectrum for the model is shown in Fig. 
The two models have log likelihood difference — logL = 127: they are easily distinguishable. 
The speed of sound has a signihcant effect on the CMB anisotropy. Figure ^ compares the 
dark matter and dark energy contributions for quintessence and fc-essence models. The small 
sound speed results in distinctive oscillations in the case of /c-essence. 

Can the effect of the sound speed be distinguished from that of other cosmological pa¬ 
rameters? There is already a large degeneracy in these parameters, so it would not have 
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FIG. 4: The CMB power spectrum for ^-essence (solid line), the model with Cg = 1 (dot-dashed 
line, — logL = 127) and the best-fit spline model (dashed, — logL = 28). Models that are dis¬ 
tinguishable can still fit quite closely. The lower diagram shows ACi/Ci (relative to the fiducial 
model) and the cosmic variance envelope (with no particular normalization). 

been too surprising if allowing a variable speed of sound merely expanded the pre-existing 
degeneracy. This problem is addressed by considering, as above, quintessence models which 
have = 1, but allowing the values of Gfe, flcDM, (quintessence) and ho to vary (subject 
to the flatness condition G;, -|- Gcdm + = !)• For the comparison models, the equation 

of state, w, is taken to be constant as a function of the scale factor a, but is allowed to 
vary from model to model. Minimizing the log-likelihood over these parameters, the best 
£t gives — logL = 32, with parameters fib = 0.05, flcoM = 0.35, flq = 0.60, ho = 0.48 and 
w = —0.78. This £t was found using well known minimization schemes. It seems likely 









FIG. 5: Comparison of the power spectra for dark matter (CDM) and dark energy (Q or K) for 
models with dark energy in the form of cosmological constant (A), quintessence {Q) and fc-essence 
{K). All models have ^cdm = 0.3, Qbh'^ = 0.02, and h = 0.65. For models with cosmological 
constant, there is no perturbed dark energy component. Note the distinctive oscillations of the 
k-essence component associated with having <C 1. 

that — log L = 32 is the best that can be done for this class of models, as the resnlt is qnite 
insensitive to the parameter valnes with which the mimization is started. Hence, the speed 
of sonnd is distingnishable. 

Thns far, onr hdncial model has been compared with two kinds of = 1 models: one 
with all other parameters, inclnding w{a), identical and one with constant w, with the 
parameters (inclnding w) adjusted so as to minimize the likelihood difference. In both cases 
the hdncial model was easily distinguishable. We next test the possibility that some other 
form of quintessence, with general w{a) and canonical kinetic energy, can reproduce the 
CMB anisotropy of the time-varying Cg model. For this purpose, we consider models with an 
equation of state given by a cubic spline. P That is, we introduce six new parameters into 
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our model: the values of w at a = 10“^, 10“^, 10“^, 10“^ and w at the two extremes of a, 
which lie at a = 10“^^ and 3.8. (Introducing more spline points has a negligible effect on 
our results.) The equation of state at other values of a is then given by a piecewise cubic 
(in logo) function whose coefficients are chosen so that it passes through the points selected 
above and has a continuous second derivative. We still allow fib, flcoM, fig and to vary 
(again enforcing flatness), and now allow, for completeness, the spectral tilt Ug to vary as 
well. The model therefore has a total of ten free parameters. The minimum log-likelihood 
difference found was — logL = 28, which is not a signihcant improvement. This model has 
fife = 0.05, flcDM = 0.34, VLq = 0.61, ho = 0.47, ng = 1.02. The equation of state is shown 
in Fig. ^ and the CMB power spectrum is shown in Fig. ^ (dotted line). 

The spline technique always produces a very smooth looking equation of state, compared 
to the rapidly varying equation of state (particularly near a = 0.01) given by the actual 
A:-essence model. To see that this does not affect the analysis - that the spline equation of 
state has sufficient freedom to closely mimic that for the hducial model - we compare two 
models, one with Cg = 1 and one with an identical speed of sound to the f;-essence model. 
The equations of state of these models are allowed to vary using the spline technique, but the 
cosmological parameters are hxed to be the same as for the hducial model. The minimum 
— logL for the Cg = 1 model is 70, whereas the log-likelihood difference for the model with 
the hducial sound speed is much less than one. Thus, the spline technique appears to do a 
very good job of modeling the relevant details of the equation of state. 

In sum, we have seen from our example that it is possible to distinguish models with 
Cs = 1 (e.g. models with canonical kinetic energy density) from models with 1, as 

occurs in models with non-canonical kinetic energy density, such as k-essence. Our studies 
show that distinction depends on c* -C 1 and Qq being at least a few percent at the last 
scattering surface so that the small value of the speed of sound ahects the acoustic peaks, 
which can be precisely measured, as well as the large-angular scale anisotropy. If Og -C 1 
at the last scattering surface, then the speed of sound only affects the large angular scale 
anisotropy and is difficult to distinguish because of the large cosmic variance at those scales. 

As a technical device for the study, we have introduced a cubic spline method for studying 
varying a general equation of state w{a) with a hnite number of htting quantities. The same 
technique can be extended to include general time-varying Cg. In this way, near-future 
observations of the CMB may be used to constrain general models of quintessence without 
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introducing strong priors into its nature. 
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